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Abstract. Let g be a fixed modular form of full level, and let {fj k} be a basis of holomorphic 
cuspidal newforms of even weight k, fixed level and fixed primitive ncbcntypus. We consider the 
Rankin-Selberg L-functions L(l/2 + it, /wj. (gig) and compute their second moment over ixT and 
k x K. For if 3 / 4 + e ^ T < K$/ A - £ we obtain an asymptotic formula with a power saving error 
term. Our result covers the second moment of L(l/2 + it + ir, fj k )L(l/2 + it — ir, fj^t) f° r am/ 
fixed real number r, hence also the fourth moment of L(l/2 + it, fj,k)- For the proof we develop 
a precise uniform approximate functional equation with explicit dependence on the archimedean 
parameters. 



1. Introduction 

Moments of families of L-functions on the critical line play an important role in arithmetic. Often 
they constitute a crucial input for subconvexity or non-vanishing results that have far-reaching 
applications. In this article we consider degree four Rankin-Selberg convolutions of a fixed modular 
form with holomorphic cusp forms of large weight k and their associated L-functions at a point 
1/2 + it high on the critical line. We aim at an asymptotic formula for the second moment with a 
power saving error term, where k and t are both large, but possibly in different ranges. A special 
case treats the fourth moment of L-functions associated with holomorphic cusp forms. 

More precisely, we shall consider the following situation. Let N 1 be an integer and \ an even 
primitive Dirichlet character modulo N, in particular N ^ 2 (mod 4). For even k ^ 2 let 

oo 

(1) fj,k(z) = J2 A J , fe (m)m i ^ i e(mz), 1 < j < 9 k (N, X ) ~ dim c S k (N, X ), 

be an orthogonal basis, consisting of arithmetically normalized newforms, for the cuspidal space 
Sk{N, x)', such a basis exists by the assumption that \ is primitive. It follows from [CH, Theoreme 1] 
that with the exception of finitely many pairs (N, k) we have (see the end of this introduction for a 
precise statement) 

(2) e k {N, X ) =kN 1+ °^. 
For a cusp form / € S k (N, \) we define 

r(fc - 1) 

P[J) - (4 7 r) fc - 1 ||/|| 2 ' 

where 

||/|| 2 :=/ \f(z)\ 2 y k - 2 dxdy 
Jr (N)\u 
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is the usual Pctcrsson norm. For the Hecke eigenform fj k it is known by the work of Iwaniec and 
Hoffstcin-Lockhart [Iwl, HL] that 

(3) p(fj, k) = (kN)- 1+ °W. 

Let g be a fixed automorphic form for the full modular group, which can be 

• a holomorphic cusp form of some weight k' < k with Hecke eigenvalues A g (n) ; 

• a Maafi cusp form of weight zero, spectral parameter r = \f\ — 1/4 ^ 0, sign e g € {±1}, 
with Hecke eigenvalues A s (n); 

• an Eisenstein series E r := E(z, 1/2 + ir) with r £ R \ {0}, or Eq := j^E(z, s)| s= i/2, with 
Hecke eigenvalues X g (n) = ~}2 ab=n {a/b) lr . 

Let 

(4) L(sJ j , k ® 9 )=L(2s, X )f; Xj ' kin) ^ n) 

* — ' n s 
n=i 

be the corresponding Rankin-Selberg L-function. Our assumptions on fj k and g imply that for 
g cuspidal this is really the usual Rankin-Selberg L-function (see [Li], in particular Example 2 on 
p. 146). This L-function is entire (and by a deep result of Ramakrishnan [Ral] it is even cuspidal 
on GL4 over Q, but we will not need this). In the case g = E r we take the right hand side of (4) 
as the definition of the left hand side, then by the multiplicative properties of Hecke eigenvalues we 
obtain 

(5) L(s, f hk ® E r ) = L(s + ir, fj,k)L(s - ir, f jik ). 

Let W\,W2 : (0, 00) — > [0, 00) be fixed smooth functions with nonempty support in [1,2], and let 
T, K ^ 1 be two sufficiently large parameters. The aim of this article is to obtain an asymptotic 
formula for 

f°° / t\ f k - 1 \ 0k{ - N ^ 

(6) I(T,K):= WA-) ]T W 2 ( — ) £ p{f jM )\L{l/2 + it,f j , k ®g)\*dt. 

J ° ^ ' fc=0(2) ^ ' J = l 

The square root of the analytic conductor of |L(l/2 + it, fj. k ® g)\ 2 is approximately given by (cf. 
(24) below) 

AT 2 / k 2 

For j G No and rGlwe shall need the following smooth averages: 

£j ( T > K y-=YK\ \ Wl r ) Wi \k) log " c(< ' x) dt dx ~ logJ {T + K) ' 
M ir (T,K) : = — J J wi (-) w 2 [-)c{t, x y r dtdx « 1. 

The shape of the asymptotic formula depends on the type of g, hence we formulate our main 
results in three separate theorems. We derive an asymptotic formula which is nontrivial precisely 
when 

(9) ^ T ^ K 5/4-e 

for some e > 0. 

Theorem 1. For g cuspidal there are constants ao,ai € R depending only on N and g such that 
X(T,K) = TK(a 1 £ 1 (T,K) + a o C o (T ) K))+O If 0,w 1 ,w a ,e {(TK) 1+e (T 4 K- 5 +T~ 3 K i )) . 
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The leading constant is given by 

ai = ±L(l,Ad 2 g)l[(l-p- 2 ). 

p\N 

Theorem 2. For g = E r with r ^ there are constants ao, 0,1,0,2 € R and b± € C depending only 
on N and r such that 

I(T, K) =Tk(^2 ajCj(T, K) + b + M ir (T, K) + 6_M_ ir (T, K) 

+ N , r , WltWa , e ((TA') 1+£ (T 4 ^- 5 + T- 3 K 4 )) . 
XVie leading constant is given by 



a 2 = i|C(l + 2ir)| 2 n(l-f- 2 ); 



p\N 



moreover, 



6± = ic(l±2ir) 4 n(l-p- 2T4ir )- 

p|JV 

Theorem 3. For g = Eq there are constants clq, 01, 02, 03, 04 € R depending only on N such that 



1(T,K) =TK\Y^a j C j {T,K)\ + N , Wl , Wa , E ({TK) 1+£ (T 4 K~ b + T~ 3 K 4 )) 
S'=o ' 
The leading constant is given by 



p\N 

Remark 1. For T — K and fixed W\ t 2 the above asymptotic formulae take a particularly simple 
shape as Cj(K, K) is a polynomial in \ogK of degree j and M.i r (K, K) is proportional to K 4lr . 

The present paper was inspired by recent work of Kim-Zhang [KZ] based on a similar result of 
Duke [Du] which in turn was motivated by a paper of Sarnak [Sal]. These three papers have in 
common that they estimate an archimedean family of L-functions of conductor T 8 and size T 2 [Sal] 
(resp. T 3 [Du, KZ]), without using approximate functional equations and orthogonality properties 
of the family. Rather, the argument is based on L 2 -techniqucs, and the idea is to express the mo- 
ment as the norm of an automorphic form on hyperbolic 3-space [Sal] (resp. 5-space [Du]) when 
restricted to a certain cone. This very elegant and flexible approach gives upper bounds with the 
correct log-power, but (in its basic form) does not yield an asymptotic formula of the above kind. 
Moreover, it requires the various archimedean parameters to be of the same size. With more work, 
these problems could be dealt with essentially by treating the off-diagonal contribution nontrivially 
as in [PS]. We also note that there are many other results in the literature on various moments of 
Rankin-Selberg L-functions, e.g. [HM, LLY, Sa2]. Particularly interesting is the preprint [Yl], where 
an upper bound for the GL3 x GL2 analogue of the quantity (6) for T w K is obtained. The paper 
[Kh] is also somewhat similar in spirit to the present work. General conjectures for a wide variety 
of moments of L-functions have been formulated in [CFKRS], sec in particular Conjecture 4.5.1, 
where the family of Theorems 2 and 3 is considered for fixed t and k. It seems very complicated to 
integrate the right hand side of [CFKRS, (4.5.8)] over t and k. 

We use essentially the same type of family as in [Du, KZ]. The novelty in our work is the derivation 
of an asymptotic formula with independent parameters T and K] there are not many results of this 
type in the literature. An upper bound of essentially the right order of magnitude for T(T, K) in 
Theorem 1 follows immediately from [JM, Section 3], even without integration over t. The main 
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point here is really to prove an asymptotic formula with a power saving error term which turns out 
to be quite strong in the case T — K. Our family has conductor (T + K) a and size TK 2 , hence 
in the extreme cases of Theorem 1 (cf. (9)) we are dealing with a family of conductor K 8 and size 
LT 11 / 4 , or a family of conductor T 8 and size T 13 / 5 . Theorem 3 states a precise asymptotic formula 
for the fourth moment of certain modular L-functions (cf. (5)-(6)), improving substantially on [Du, 
Theorem 1] and [KZ, Theorem 2]. 

Our approach is based on classical techniques, in particular an approximate functional equation, 
Petersson's formula, Voronoi summation etc. It may be observed, however, that at least in the 
range (9) we do not have to treat shifted convolutions sums explicitly, although of course they 
appear throughout the discussion of the off-diagonal term. In addition, we do not have to use any 
nontrivial bound for twisted Kloosterman sums. 

As an interesting feature we remark that in the Eiscnstein case the off-diagonal term contributes 
to the main terms in Theorems 2 and 3, and it requires a non-trivial manipulation to recover the 
exact shape of this off-diagonal contribution. The shape of the asymptotic formula in Theorem 2 
is particularly interesting, as it contains oscillating secondary terms A4 ±i r (T, K). This seems to 
be a new phenomenon of this family. The fact that the off-diagonal term contributes towards the 
main term seems to be a general feature in moment computations of Rankin-Selberg L-functions 
with Eisenstein series, see for instance [Bl] or [DFI]. It would be desirable to have a conceptual 
understanding of this empirical phenomenon. It also indicates that there is some intrinsic difficulty 
in deriving an asymptotic formula for a fourth moment of automorphic L-functions as in Theorems 2 
and 3. 

Our method could in principle treat more general assumptions on the parameters of fj t k and g, but 
this would infer technicalities that may obscure the paper rather than be useful in applications. It 
is possible by standard techniques to remove the harmonic weight p(fj,k), see e.g. [IK, Section 26.5]. 
It should also be possible to replace the family of holomorphic cusp forms of weights k x K by a 
basis of Maafi forms of spectral parameters tj X K. In this situation, however, the analysis would be 
somewhat different, and in certain ranges one would encounter effects where the conductor drops, 
see e.g. [Y2] for a related discussion. Finally and most importantly, it is possible to prove asymptotic 
formulae in the spirit of Theorem 1 in other ranges, for example in K E $J T $J K , and most likely 
also for T considerably bigger than K 5 ^ 4 . This would require, however, a careful treatment of shifted 
convolution sums, in the latter case of large T with oscillating weight functions. As remarked above, 
the present approach avoids shifted convolution sums. 

In the course of the proof we use two auxiliary results that may be useful in other situations. 
These are certainly known to specialists, but we could not find them in the literature. Proposition 1 
proves a relatively sophisticated form of the approximate functional equation, while Proposition 2 
establishes a Voronoi summation formula for Eiscnstein series L r . 

Finally we remark that [CH, Theoreme 1] can be used to show that 8k(N, x) = for at least one 
even primitive x if an d onr y if the pair (N,k) is one of (1,2), (1,4), (1,6), (1,8), (1,10), (1,14), 
(5,2), (7,2), (8,2), (9,2), (11,2), (12,2), (13,2), (15,2), (17,2), (19,2), (21,2). In other words, (2) 
holds for (TV, k) outside this explicit set. 

Notation and conventions. In an effort to lighten the notational burden, we will use the following 
conventions. The symbol e denotes an arbitrarily small positive constant and the letter A denotes an 
arbitrarily large positive constant, not necessarily the same at each occurrence. All implied constants 
may depend on e, A, Wi, W<i. This allows us to write, for instance, C e logC <C C e . All implied 
constants may also depend on N and the parameters of the fixed modular form g; the reader may 
check that the dependence on these quantities is always polynomial. The phrase "negligible error" 
means an error of size 0((T + K)- A ). By (9) this is the same as 0(T A ) or 0(K A ). 

Acknowledgements. We thank the Institute for Advanced Study for inviting us to the inspiring 
"Workshop on Analytic Number Theory" in March 2010, where good progress was made on this 
paper. We also thank the referee for a careful reading and numerous comments. 
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2. Preparatory material 

2.1. A uniform approximate functional equation. We start with a general result that may be 
of independent interest. In order to prove an asymptotic formula for X(T,K), it is very important 
to have an approximate functional equation with a weight function that is independent of the 
archimedean parameters. An approximate functional equation of this kind was developed in general 
in [Ha, Theorem 2.5]. Here wc need a slightly stronger variant of this result with a better error term, 
which forces us to use a somewhat more complicated main term, similarly as in [HB, Lemma 1]. We 
shall prove a general result since it requires the same amount of work as the proof of the special case 
of Rankin-Selberg L-functions we are interested in here. 

At this point we take the opportunity to correct an inaccuracy in our earlier work [BH]. A uni- 
form approximate functional equation as stated in Proposition 1 below should have been used in 
[BH, (2.12)]. This would justify, at the cost of an admissible error O e , A (D l / 2+e T- A ) , the tacit 
assumption that the weight function V in [BH, (7.2)] is essentially independent of t. 



We keep the notation of [Ha]. Let F be a number field of degree d, and let n = ® v tt v be an 
isobaric automorphic representation of GL m over F (cf. [Ra2]) with unitary central character and 
contragradient representation tt. The corresponding L-functions are defined for 3?s > 1 by absolutely 
convergent Dirichlet series as 



oo oo 



n—l n—1 

which extend to meromorphic functions on C with finitely many poles, and these are connected by 
a functional equation of the form 

N s/2 L{s, tt 00 )L(s, tt) = kN {1 ~ s)/2 L{1 - s, froo)L(l - s, f ). 

Here Af is the conductor (a positive integer), k is the root number (of modulus 1) and 

3 = 1 \ ^ / j=1 

for certain fj,j G C which satisfy 
(10) Ztfij > 



to 2 + 1 2 

by a result of Luo-Rudnick-Sarnak [LRS] . We put 

1 . N -II \ £(l/2,7Too) 

r?i = — -, ri := mm n,- , A = — ; — ; h 

13 4 2' ' i«w" Jh L(l/2,7Too) 

and we define the analytic conductor (at s = 1/2) as 



N j-r 1 N j-f . 



For a multi-index n S Nq wc write |n| := n(l) + ■ • • + n(2md) and 

(11) f/- n := n»?7 n(JW -%- n(!W) . 

By a result of Moltcni [Mol, Theorem 4], the coefficients of L(s,ir) satisfy the uniform bound 

(12) ^ |a n | < e z 1+£ C £ . 

niCx 

To be precise, axiom (A4) in [Mol] includes ^ftfij ^ but (10) is sufficient for the proof. 
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Let us assume that L(s,ir) is entire, then we can state the following result: 

Proposition 1. Let Go ■ (0, oo) — > M be a smooth function with functional equation Gq(x) + 
Gq{1/x) = 1 and derivatives decaying faster than any negative power of x as x — > oo. Let M £ N. 
There are explicitly computable rational constants c n j € Q depending only on n, £, M , m, d such 
that the following holds for 

(13) G(x) := G (x) + ]T c^t?"" (x— ) G (x). 



0<\n\<M 
0<£<\n\+M 



a) The function G is smooth and its Mellin transform G is holomorphic everywhere except for a 
simple pole at s = with formal Laurent expansion 

(14) G(s) ~ - + J2c jS \ cj = J a G'(x)(\ogxy +1 dx = d j + Oir 1 )- 
Moreover, one has 

(15) -^G(x) < (l + x)- A ; G{s) < (1 + \s\)- A for \Us\ < a and \s\ > 1. 

ox 1 

Here dj and the implied constants depend at most on j, A, o~o> M , m, d, and the function Go- 

b) For any e > one has 



(16) L(1 M = £ ^G ( » + kA £ ^ G ( » + O^G^) 



n— 1 v \ v / n— 1 v x v 

where the implied constant depends at most on e, M , m, d, and the function Go- 

Remark 2. For sclf-contragradicnt representations 7r we can strengthen Proposition 1 by imposing 
the additional symmetry Cn,e = Cn,t for the involution n H> n defined by rj~ n = rj~ n , that is 
n(2j — 1) := n(2j) and n(2j) := n(2j — 1). Then G defined by (13) is real-valued and the two 
main terms in (16) are identical. Indeed, in the self-contragradient situation k = X = 1 and the 
coefficients a n are real, hence we can replace the original G by its real part without affecting the 
validity of (14)-(16). This corresponds to replacing the original coefficients c n .i by [c n .i + Cn l e)/2, 
upon noting that the original coefficients arc rational. 

Proof. Wc follow closely the proof of [Ha, Theorem 2.5]. We can write Go as an inverse Mellin 
transform 

(17) G (x) = ± f x- s H(s)-, v>0, 
where H (s) := sGo(s) extends to an entire functions satisfying 

(18) H(0) = l, H(s) = H(-s) = H(s), H{s) < (1 + \s\)~ A for < a . 
As in the Erratum of [Ha] wc define 

L (S, TToo) .= -G Al — — — \TlMO T + G ' 

then by the proof of [Ha, Theorem 2.1] we have the exact formula 



L(i/2, tt) = V ±w ( + aV ( ") 

vv^y l^j 
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with the weight function 

W{x) = ^-( x-C-^Ffair^Hia) — , o > 0. 

2?" J ( CT ) s 

Now the idea is to approximate W by a function G of the form (cf. (11)) 

G( x ) = ^-f x- s (l+ ^-"(-^W)-, a>0, 

71-4 ^(°") ^ 0<\n\<M ' 

0<t<\n\+M 

which, by (17), is precisely the function defined in (13). In fact 

(19) G(s) = (l + c n jJi- n {-') t )&o{'), 

^ 0<|n|<M ' 

0<£<|n| + A/ 

hence (14) and (15) are immediate from (18) and Gq(s) = H(s)/s. The rapid decay of W together 
with (12) implies (16) as soon as the constants c n ,£ are chosen such that 

W(x) = G{x) + M , m .Av~ M )- 

We shall derive this from 

W{x) - G{x) = / x- s (c- s ' 2 F{ Sl ttoo) - 1 - Y, c r^V- n (sy) H{s) — 

m ^ 0<|n|<M ' S 

0<e<\n\+M 

by showing that for explicitly computable rational constants c n ^ £ Q we have 

(20) C-^ 2 F(it,n 00 ) = l+ Y, ^ eV - n (-ztY + M ^. d (v~ M {\t\ + \t\ 2M )), t £ M. 

0<|n|<Af 
0<£<|n|+A/ 

Note that this strengthens [Ha, Lemma 4.1]. 

In proving (20) we can restrict to the range \t\ < rj, because for \t\ r\ the approximation is 
trivial for any constants c n j. Indeed, rj 3> m 1 by (10), hence we have 

1 «M. m \t\ M < r,~ M \t\ 2M , 
and similarly for < \n\ < M and < £ < \n\ + M we have 

cn^i-it) 1 v~ ]nl \t\ £ «M, m 77" |Tl| |i| |n|+M «S r7- M |t| 2M . 

Let us now assume that \t\ < 77, then our starting point is the identity 



C-^F(it,7T 00 ) = Ull 



- u r(i + f + |)r(i + f ) 



2 2 



1 1 Mj 
4 2 



1 1 

2 + 2 



r(i + f -f)r(i + f) 

"id /-t/2 /tV 



exp \^Y J Q (~f + * T ) ~ lo S '^') dT 



In the integral we have |ir| < 77/2 and \r]j + ir| > 77/2, hence the well-known asymptotic expansion 
^(z) = logz- Y ^ + 0*,m(\z\- M ), SOOO, 

0<n<A/ 

where B n £ Q is the n-th Bernoulli number (cf. [GR, 8.361.8]), yields 

/ t/2 (y (Vj + ir) - log%) dr= Y + °M (|i| + |i| M+1 )) 

\ ' 0<n<M 



s 
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for some constants b n ^ £ Q. It follows that 

md „t/2 /p, \ 



i d 



J2 Km" E for + M)^") + o M>m4 (r M (1*1 + i*i m+1 )) 



0<n<M j=l 

Finally we approximate uniformly the exponential function for imaginary arguments by its Taylor 
polynomial of degree M — 1 to arrive at (20) for some c n j S Q depending only on n, ^, M , m, d. 
Here we use that for M \n\ ^ M 2 and 1 ^ £ ^ |n| + M we have 

rT W |^ ^ rr H (|t| + |t| H+M ) «A/, m V~ M (1*1 + \t\ 2M ) , 

because r\ ^$> m 1 by (10) and |i| < 77. □ 

We shall fix once and for all a weight function Go : (0, 00) — > R as in Proposition 1, so we will not 
display the dependence of our statements on this function. We shall fixJ/eff later in the paper. 
For t real and k 2 an integer we shall apply Proposition 1 to the particular L-function 



£t,j,fc(s) := L(s + it, fj. k ® g)L(s + it, f jtk (8 

This L-function is associated with a self-contragradient representation tt, namely the isobaric sum 
°f fj,k ®g <8> I det 1 14 and its contragradient. According to Remark 2 we can and we shall assume that 
G is real-valued, then (16) becomes 

(21) \L(l/2 + it, f hk ® , 9 )| 2 = 2 £ ^=G \^=j + O^C 1 ^) . 

We now compute all the necessary data occurring in this situation. We have 

(22) m = 8, d = 1, TV = TV 4 . 
Using (4), we find that the Dirichlet coefficients of £t,j,k(s) are given by 

(23) a n = E x(di)x(d2)d^ 2i Uf t X hk {m 1 )X j:k (m2)Xg(m 1 )X g {m 2 )m^ t mf, 

so that (12) is satisfied by previous remarks or by standard bounds. There are constants V\, . . . , 1/4 € 
C depending only on g such that the archimedean parameters fxi, . . . , fxs € C of £t,j k(s) are given 
by the numbers + Vj + it and their complex conjugates. Hence we have 

k $lvj i(^Vj +t) , 



(24) ^4 4 



1 



(k+\ t \f 



2.2. Voronoi summation. The purpose of this section is to compile summation formulae for the 
Hecke eigenvalues X g (n) twisted by a finite order additive character. This generalizes Voronoi's 
original formula for the divisor function (without twist). 

Proposition 2. Let a and c be coprime positive integers, and let F : (0, 00) — > C be a smooth 
function of compact support. Then 

(25) cf; X g (n)e (n^j F{n) = X g („) £ e (rnfj J™ F(x) jf dx , 
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where 

J+(x):=2iri k J k - 1 {x), J g (x):=0 
if g is a holomorphic cusp form of level 1 and weight k; 

(26) J+(x) := ~J j Y 2ir {x) + Y_2ir(x)), J~(x) := e 9 4cosh(^r)i^ 2jr (x) 

M cosh(7rr) ' y 

if g is a Maafi cusp form of level 1, weight ; Laplacian eigenvalue 1/4 + r 2 , and sign e g G {±1}- 
For g — E r (r£8) the same formula holds with as in the Maafi case (with e g = 1), except that 
on the right-hand side the following polar term has to be added: 



. , r°° , x s ±ir 

X)C(l±2<r) J o (J) F(x)dx 



(27) > C(l±2ir) / l-r) F(x)dx forr^O, 



x 



(28) J [log ^) + 2 7 J F(x) dx for r = 0. 

Proof. For Maafi cusp forms the result is due to Meurman [Me, Theorem 2]; for holomorphic cusp 
forms it is due to Jutila [Ju2, Section 1.9] and Dukc -Iwanicc [DI, Theorem 4]; for g = Eq it is due 
to Jutila [Jul, Theorem 4]; for g = E r (r ^ 0) we provide the proof below. 
Let g = E r (r^O), then the Hecke eigenvalues can be explicitly described as 



ab—n 

We shall derive the Voronoi formula from analytic properties of the Dirichlet series 



D(g, x, s) := \ g (n)e{nx)n s , 5fts > 1, 

71=1 

summarized by the following lemma. 

Lemma 1. For (a,c) = 1 the Dirichlet series D (g, 2 ,s) can be analytically continued to a mero- 
morphic function which is holomorphic in the whole complex plane up to simple poles at s = 1 ±ir, 
where the residues equal £(1 ± 2ir)/c 1±2lr . It satisfies the functional equation 



(29) D [g,-, sj =tt- 1 2 Zs - l ^-J r(l - s + ir)T(l - s - ir) 



l-2s 



x < cos(irir)D I g, — , 1 — s I — cos(7rs)D g, , 1 



c 



Remark 3. The limiting case r — s- of this lemma is due to Estermann [Es] and served for Jutila 
as the starting point for his Voronoi formula, see [Jul, Lemma 1]. Correspondingly, the limit of (27) 
under r — >• equals (28). 

Proof of Lemma 1. The result follows from [Mot, Lemma 3.7] upon noting that the Dirichlet series 
D (g, §, s) is a shift of a so-called Estermann zeta-function: 

oo 

D (.9, -,s-irj = o- 2lr {n)e (— ) 

n—l 

The proof of [Mot, Lemma 3.7] is based on analytic properties of the Hurwitz zeta-function. In the 
Appendix we provide an alternate proof of Lemma 1, based on the modularity of g = E r , which 
displays the similarity with the cuspidal case. □ 

With Lemma 1 at hand it is straightforward to deduce (25) with the additional polar term (27) 
on the right hand side. The left hand side of (25) equals 

cjr\ g (n)e(n-)F(n) = ^ f F(s) D (g,-,s) ds, 
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where F denotes the Mellin transform of F. We shift the contour to the vertical line at —1 and 
record the contribution of the residues at s = 1 ± ir; we obtain 

£ X 9 (n)e(n^F(n) = £ ± ir) + £ / *(«) ±s) *■ 



c 

n=l 



We observe that the sum on the right hand side equals (27), hence we are left with proving that the 
remaining term equals the right hand side of (25). By (29) the term in question can be rewritten as 



If f W (fpr(i 



s + ?r)r(l — s — ir) 



'(-i) 

x |cos(7rir)Z) (a, — , 1 — s \ — cos(ns)D (g, — — , 1 — sj | ds. 

We apply the change of variable s — > 1— | and unfold the series D (a, ±-, s). By absolute convergence 
we see that the previous display equals 

, s ( o\ 1 /" f2n;Jri\~ s , s /s + 2ir\ /s-2ir\w s\ 

e a » ( + - c) 2^ y (2) (-r-j r (— ) r (—j f i 1 - 2) ds 



n=l 



a\ 1 /" { 2ny/n\ s /7rs\ /s + 2ir\ /s — 2ir\ ~ / s 



By [Er, 6.8.17 & 6.8.26] and (26) this is the same as 

tv-)E.H)^/ M (^)^?C)f(i-3)*. 

which is precisely the right hand side of (25), thanks to the following identity for t > 0: 

JL rv>) P (1 - 1) * = ^ ^ r . FWl -/ 2 dI } „ 

/•<» 

= / Jf (tVx) dx. 

Jo 

The proof of Proposition 2 is complete. □ 

2.3. Bessel functions. We start with some standard bounds; more precise results can be found in 
[HM, Appendix]. It follows from the power series expansion that 

k— 1 

(30) J k -i(x) « ^ttt , 0<x^l, 

T{k) 

uniformly in k £ N. It follows from the asymptotic formula and the power series expansion that 

(31) Jk-i(x), Y 2lr {x) < x- 1/2 , x > 0, 
for fixed k £ N and rgR, Finally we recall the bound 

(32) K 2lr (x) « e- x , x>l, 

for fixed r £ R, which follows again from the asymptotic formula. 
For a smooth function h : (0, 00) — > C of compact support let 

h(t):= f h(x)e(xt)dx, hit) := — L / hUfx)e ( — ) -^L 
Jo V27T Jo V 27r / V x 

Then the following summation formula holds. 
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Lemma 2. Let K > 0, £ > 0, and h : (0,oo) — > R a smooth function of compact support. Then 

E <-* ( V) J '-«> ^ ( e (i " i) 4 (f ) ) + ° (w £ ^ * 



fc=0 (2) 

ura'i/i an absolute implied constant. 

Remark 4. For odd weights one has the similar formula 



fe=l(2) 

Remark 5. It is well-known that an individual Besscl function Jk~x(x) is hard to control in the 
range k <C x <C k 2 . The lemma shows, however, that on average over k its values cancel almost 
completely until x ~ fc 2 , when the asymptotic behavior becomes stable. For further discussion of 
the error term see [Iw2, p. 87]. 

Proof. This is based on [Iw2, Lemma 5.8] and can be found in [Kh, Lemma 2.3]. □ 

Finally we state a result which is useful for estimating the Hankel-type transform occurring in 
the Voronoi formula (25). 

Lemma 3. Let F : (0,oo) — > C be a smooth function of compact support. For s G C let B s denote 
either of the Bessel functions J s , Y s or K s . Then for a > and j € N we have 

J F(x)B s (ay/x)dx = ± l—J J ^— (F(x)x~i)x^ J ~ B s+j (ay/x) dx. 

Proof. The Bessel functions B s satisfy the recurrence relation (x s B s (x))' = ±x s B s _i which trans- 
lates to 

(aVx~yB s (aVx~) = ± ^^-((aV^) s+1 B s+1 (ay/^)). 
Using this identity and applying integration by parts j times we obtain 

J F(x)B.(ay/x)dx = ±\-p\ J F(x)(aVxy s -g-{{aV^) s+3 B s+J (a^))dx 
= ±(^J J —(F(x)(aVx~)- s )(aVx~y + iB s+] (aVx~)dx 

( 2 \ j f°° 8 j a a + i 

= ±(-j / ■g—(F(x)x-2)x—B s+: j(ay/x)dx. 

□ 

2.4. Fourier coefficients of cusp forms. For a Dirichlet character x modulo N and a positive 
integer c divisible by TV let 

b x {m,n,c) := ^ XW e 

d (mod c) 

be the twisted Kloosterman sum. We formulate Petersson's formula for the basis of holomorphic 
cusp forms specified in (1). 

Lemma 4. For k, m,neN and k 3 we have 

e k (N. x ) a , , , 



J2 p(/Afc)A A fc(m)A J)fc (n) = 5 m ,n + 2™~ fc £ xl ' } J k - 

3=1 W|e 

Proof. This follows from [IK, Proposition 14.5]. □ 
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Wc shall often use the following standard result. 
Lemma 5. We have the uniform bound 

(33) J2 IV™)| 2 M 1+e - 

Proof. See [IK, (1.80) & (14.56)] for stronger bounds. □ 

3. The diagonal term 

We shall assume (9) for the rest of the paper, since otherwise the asymptotic formula is trivial. 
We substitute the approximate functional equation (21) in the special case (22)-(24) into (6) and 
choose M := 3. We note that both G = Gt,k and C = Ct : k depend (mildly) on t and k, but in the 
support of Wi and W 2 we have C t ,k ~ C := (T + K) 8 . By (2)-(3) (or directly by Lemma 4) the 
error term contributes at most 

(34) < TKC 1/4+E (T + K)- M < (TK) 1+E (T + K)~ l . 
After inserting the main term, the j-sum in (6) equals 

( 35 ) 2 E G ^J^ E ^r^u^^—lT E p(/,, fc )A j ,,(m 1 )A i , fc (m 2 ). 
By Lemma 4, the innermost sum equals 

S x (mi,m 2 ,c) f 4ir y /m 1 m 2 



(36) Wh+^E ^ 

AT|c 

The (5-term contributes 



Let G denote the Mellin transform of G and similarly for other functions. By [Li, p. 145] we have 

|A g (m)| 2 L{s,g(g>g) 



E 



C(2s) 



when 5 is cuspidal. If g = E r is an Eisenstein series, we take the preceding display as the definition 
for L(s,g®g), so that L(s, E r <E)E r ) = (^(s) 2 £(s + 2ir)((s — 2ir). With this notation the contribution 
of the 5-term in Petersson's formula equals 



J ° V 7 fc=0(2) V 7 

x-L /" £(1 + 2»t + 2s, X )L(1 lit + 2.s, x) ^t^tj f ^ g s/2 G( S ) d S dt. 



(37) 

x 

2m j (1 

We shift the contour to the line SRs = — 1/4 + e, then the new integral contributes (cf. (15)) 

(38) < (TK) 1+£ {T + A') -1 . 

There is pole at s = whose order u depends on g: if g is a cusp form then v = 2, if g = E r 
with r / then u = 3, while for g = Eq we have v = 5. The residue of the pole is given by a linear 
combination (with coefficients depending at most on N and g) of 

r W 4^) E W 2 (^\L^\l + 2zt lX )L^\l~2it,x)\oeC{t,k)dt 

(39) Jo \ T J k ^ {2) \ K J 

+ 0((TK) 1+£ (T + K)- 1 ) 
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for ji + ji + j 3 ^ v — 1. The error term comes from approximating the square root of C = C tt k m 
(24) by (7) and inserting (14). The constant in front of the leading term {31,32,33) = (0,0, v — 1) 
equals (cf. (14)) 

! £ ^C(2) g) ' 9 cuspidal, 
^-^Sr 1 , 9 = E r , r^O, 
192C(2) ' S = E . 

For x trivial (i.e. N=l) there are also poles at s = ±it, but their contribution is (cf. (15)) 
(41) < (TA') 1+e T~' 4 . 

If g = E r with r 7^ 0, there are two additional simple poles in (37) at ±ir with residue 

'k-V 



(42) 



2 rw~) e ^ 

J ° W fe0(2) 



A" 



A(l + 2g ± 2»r, x)A(l - 2g ± 2jr, X )C(1 ± 2ir)\{\ ± 4ir) c ± lr /2£ (± ■ )d , 
2((2± Air) 1 lT) ' 



Here again we can approximate C 1 / 2 by (7) and G(±ir) by Go(±ir) at the cost of an error term 
0((TK) 1+£ (T + AT)" 1 ), cf. (24) and (19). 

We continue with the analysis of (39). Applying Poisson summation it is straightforward to sec 
that the main term in (39) equals, for any A > 0, 



(43) 2 



Jo 



Wi ( - J W 2 (-J A^(l + 2it, x)A fe) (l - 2i*,5c) log j3 C(i,x) d£ da: 
+ 0((TA) 1+£ A- A ) . 
Now we use the simple approximation 



which follows in a standard fashion from van der Corput's lemma (see Lemma 4.10 and Theorem 4.11 
in [Ti]). We substitute this into (43) with M := e 6 NT (note that j u j 2 < 4 and \t\ ^ 2T). Then 
the f-integral equals, uniformly for x x A, 



V xMx(n)(-iog^(-iog^ r- flWc^^H^cft + oarJO-). 

~! ran J n \7/ V 771 / 

The diagonal and error terms together contribute (cf. (8)) 

(44) TK£j s (T, A) + 0(A(TA) e ) , 

where (AT) indicates the removal of the Euler factors at primes dividing N. Let us now estimate the 
off-diagonal term. Partial integration shows, for any m ^ n and x X A, 



V T / Vm/ |log(n/m)| 
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hence the off-diagonal term contributes at most 

1 



< K(TK) E J2 



mn\\og{n/m)\ 

m,n^e NT 

< K{TK f Yl 



n minfm, n — m) 

m<n^e e NT v ' ' 



(45) «K(TKY £ -L ? «K(TKY. 

m' ,n'^e 6 NT 

The same argument shows that (42) equals, up to an error term already present in (39), 

(A*\ TKKA IT ^ C W (2±4zr)C(l±2 t r) 2 C(l±4 ? r) ^ 

(46) TKM±i r (T, K) 2((2 ± ^ G (±tr). 

Remark 6. We will see in Section 4.6 that the term (46) will be cancelled by some portion of the 
off-diagonal term. We note that we could have avoided the computation of (42) and (46) in the case 
g = E r with r ^ by choosing Go as in (17)-(18) with the additional feature that H{s) is divisible 
by s 2 + r 2 . Then Go(±ir) = 0, hence by (19) also G(±ir) = 0, and the residues vanish. This trick 
was probably used for the first time in [BHM]. 

The various error terms in (34), (38)-(39), (41), (43)-(45) are admissible in Theorems 1-3. The 
main term in (44) furnishes the main term of Theorem 1 and part of the main terms in Theorems 2- 
3. We note that Section 4.6 discusses an additional contribution to the main terms of Theorems 2-3. 
The leading constant in the various cases of g follows from (40) and (44). 

4. The off-diagonal term 

4.1. Averaging over k and t. We return to (35) which is the j'-sum in (6), and substitute now 
the Kloosterman term in (36) for the innermost sum. This gives a total contribution of 



(47) 



fc=0 (2) 

X(^i) x(d 2 ) A g (mi) A g (m 2 ) \ - S x (m 1 ,m 2 ,c) ^ / ATTy/m^m£ 



ifd|mim 2 =n 1 2 N\c 



Recall that both G = Gt,k and C = Ct.k, defined in (13) and (24), depend (mildly) on t and k, 
but we have Ct,k X C := (T + K) s . Moreover, G is real- valued as we have assumed according to 
Remark 2. First we observe that (47) is absolutely convergent, due to the rapid decay of the Bcssel 
J-function near for large k, see (30). More precisely, we can truncate the multiple sum at 



n 



at the cost of a negligible error. 
Next we write 




Wi [ £ J W 2 (^-J- ) G tl k I —£= I =: 



where fix : (0, oo) 3 — > M is smooth function depending on K, T, N and g with the following properties: 
it is compactly supported in first two variables and rapidly decaying in the third variable, uniformly 
in K and T. Let 



(48) n 1 (x,y,z):=W 1 (x)W 2 (y)G 



(2n) 4 zVC 



N 2 ({xT) 2 + {\yKff , 
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Then by (11), (13), (24) the partial derivatives satisfy, uniformly in x 6 (0,oo) 3 , 
x k (p?\x) - O^(x)) < jtk (T + K)-\ j,k€ Nq. 
In order to sum over k with the help of Lemma 2, we define 

/•OO 

Q 2 (x,u,z) := / Ox^^^e^dy. 



Then fl 2 : (0, oo) x E x (0, oo) — > C is compactly supported in the first variable. Integration by parts 
shows that it is rapidly decaying in the other two variables, and moreover for 

/>oo 

(49) n 2 (x,u,z) := / n 1 {x,y,z)e i ^dy 



we have, uniformly in x € (0, oo) x I x (0, oo), 



x k (n { 2 j \x) - fl$\x)) «j, fe (T + A')" 1 , j,k e 



Applying now Lemma 2 with £ := 4n y/rfTpm^ / c to (47), we obtain the main term 
K ^l±i j°° ( l ±R2c d\dlm 1 m 2 \ 



2 In ^— ' \ T' 87T,/mi 7719 



X(rfi) x(^2) A g (mi) \ g {m 2 ) ( 1^mxm 2 \ S x (m 1 ,m 2 ,c) 
X di +2it 4" 2it m? /4+it m^ /4 - ite \ c J cV2 
while the error term in Lemma 2 infers a total error of 

Sif T i, , n, / m I^Y( m l> m 2)C)| 

<<C K4 X, |A 9 (mi)A 3 (m 2 )| g2 
(50) ^1° 

X I^MI 2 « (^T r(r4 // 4) « {TK) 1+e (g + g) 

Here it is enough to use the trivial bound |,S x (mi,TO2, c)\ < c and (33). 
In order to integrate over t, we define 



fls(w,u, z) := / Q 2 (x, u, z)e lwx dx. 
Jo 

Then £1$ : R 2 x (0, oo) — > C is rapidly decaying in all three variables. Putting 

/>oo 

(51) ft a (w,u,z):= n 2 {x,u,z)e iwx dx, 

Jo 

we see as above, that uniformly in i£l 2 x (0, oo), 

(52) x k {p!i\x) - n ( i\x)) < iife {T + K)~\ j,ke K 
In addition, we factor out S = (di,d 2 ) to arrive at 

(di,d 2 )=l 



x(di)x(d 2 ) \ g (mi) \ g (m 2 ) ( 2^/m 1 m 2 \ S x (m\,m 2 ,c) 

X ~~ : " —J-. r-r- — C lb- 
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It is convenient to introduce dyadic decompositions. Let uj : (0, oo) — > [0, oo) be a smooth function 
supported on [1/2,2] such that 

oo 

(53) ^w(x/2 J ') = 1, x^l. 

3=0 

Then we can recast the preceding expression as 



2 ^— ' ^—^ ^— ' ^— ' \ rf?771l ' 87T,/mi 7719 ' 

± Afi,Afa,G>l (*,JV)=1 TOi,m 2 AT|c \ 1 v 

(54) powers of two (di,d 2 ) = l 



X W 



mi\ / m2 \ / c \ x(^i)x(^2) A g (mi) A g (m 2 ) / 2^/mim 2 \ 5 x (mi,m 2 ,e) 

mJ^V^J^IgJ £ 2 dl<fe m 3/4 m 3/4 c ^ c l/2 



By the rapid decay of Q3 in all three variables we can restrict the summation to 

(55) Ml M 2 <^, ^ . 

Moreover, we can assume 

(56) d\m 2 = d\mi (l + 0(T £ " 1 )) . 
We will keep this in mind for later and in particular often use 

(57) d 2 2 M 2 x d?Mi. 

4.2. Interlude on character sums. For later purposes we need to transform the expression 

2 v / m7m 2 



S x (mi,m 2 ,c)e I ± 

in (54). We proceed similarly as in [Yl, Section 6]. We start by writing 

^ / 2^mim 2 \ _ ^ f_^ d\m l + d 2 m 2 \ ^ / (i/tifmi - v /^to 2 ) 2 



d\d 2 c I \ d\d 2 c 



We infer 



dfmi + d\m 2 

b x (m 1 ,m 2 ,c)e ± — 

d\d 2 c 



^( C ) Q (A A A A A A \ ( ^1^1 + d\m 2 

b x \d\d 2 rn\ 1 d\d 2 m 2l d\d 2 c)e I ± 



4>{d\d 2 c) ' ' \ d\d 2 c 

<j>(c) t-^* ( d 2 d±d 1 d 1 d±d 2 \ 

58 = — - > x{d)e[mi — Vm 2 — 

0(did 2 c) f-f V «2C die J 



d {d\d2c) 



We write 



(d±d ± d 2l c) = g, c = gh, did±d 2 =gf. 

Note that (gf,did 2 ) — 1 by (di,d 2 ) = 1. Let us fix a decomposition c — gh with (g,did 2 ) = 1. 
It is straightforward to check, using again (d\,d 2 ) = 1, that as d runs through (Z/did 2 cZ) x with 
(c?ig?±g? 2 , c) = <?, then / runs through (Z/d 2 d 2 /iZ) x with (gf z fd 2 , d 2 d 2 h) = d\, and this is a bijection 
between the relevant residue classes. With this notation we have 



d 2 d ± d\ = ±gf(gf =p d 2 )/d\ (mod d\d 2 c). 
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Hence we can recast (58) as 



<m ST ^* ( 9fT<k \ I f(gfTd 2 )/d 1 f 

w^cj f ? x {~^ r[ ±mi — dVh — 

(g,d 1 d 2 ) = l (gf Td2> dld2h)=di 

Using the change of variable / <— > f, we can summarize the preceding discussion as 



2^ mi m 2 \ ( Ud\mi - ^ d 2 m 2 ) 2 
S x (m ll m 2 ,c)e [ ±— — ) = e ( T 



c 



d\d 2 c J 4>(did 2 c) 



(59) v v- v-* { gfT<k \ (, (gTd 2 f)/d 1 f " 



gh=c f(dld 2 h) 
( 9 ,dxd 2 )=l ( g f Td2 , d 2 d2h) = dl 

For later purposes, we prove the following two lemmata. 

Lemma 6. Let di, d 2 , g, h, mi, m 2 be positive integers such that gh is divisible by N and di,d 2 ,g 
are pairwise coprime. Then 



f(dld 2 h) 
(gTd 2 f,dl d 2h)=d 1 



^ Ndid 2 r(h) 2 (h 7 (mi — m 2 )did 2 ). 



Proof. We use first that % is a character modulo N which divides gh. Writing 

N 

h = h\h 2 , h\ 



we see that the character sum in question equals 

( gfi T d 2 \ v^* ( t \ f 

\ x \-^r) ^ e v (TOl " W2) ^ 

h(d\d 2 hi) f(djd 2 h) 
(g^d 2 f l ,d\d 2 h 1 )=d 1 {g^d 2 f,d\d 2 h)=d 1 

f=h (d?d 2 fei) 



a j u . ._ 9Td 2 h 



With the notation 

/ =: fx + kdfd 2 hi, g\ : 
this becomes 



fi(dfd 2 h 1 ) k(h 2 ) 
(g^d 2 f 1 ,djd 2 h 1 )=d 1 (/i + fcd?d 2 ?ii,/i2) = l 

(gi=Ffcdid2?ii,?i2) = l 



The fc-sum equals 



V" n(h)lJ,(t 2 ) V* e ( (mi - m 2 )did 2 -^- ) 



il\h 2 k(h 2 
h\h 3 kd\d 2 h l = -f l {i 1 ) 

kdid\hx = ±g\ (£2) 

It is straightforward to evaluate the inner sum explicitly, but for our purposes it suffices to record 
that it is at most (h 2 , (mi — m 2 )did 2 ) which gives the lemma. □ 

Lemma 7. For r€l and Sis > 1/2 let 



(<5,AT)=1 ^ z; gfi=0 (TV) y ^ Zy y f{d\d 2 h) V 

(<*l,d 2 )=l (ff,did 2 )=l {gfTd 2 ,d 2 1 d 2 h)=d 1 
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where x is an even primitive Dirichlet character modulo N and 

X:=2 + 2s + 2ir, fx := 2s - 2ir. 

Then 

(61) S±(s) = C(2s - 2ir)C(l + 2s + 2ir)N~ 2s+2vr - p- 2 ~ 4lr ). 

p\N 

Proof. By the remarks below (58) the innermost sum in (60) can be rewritten as 

E* <^)= E* m 

f (d\d 2 h) x d(d 1 d 2 gh) 

(gfTd 2 ,d 2 1 d 2 h)=d 1 (d x d±d 2 ,gh)=g 

On the right hand side the condition on d only depends on d mod gh, therefore 

ton x(^)= E* x(d). 



{gfTd 2 ,d'id 2 h)=d 1 (d 1 d±d 2 ,gh)=g 



With the notation 



we can now write 



XW 

d(gh) 
(d 1 d±d 2 ,gh)=g 



s *(*) = E E jd^y: E ^tx^h^M- 

(&,N) = 1 (did 2 ,N)=l v 1 gh=0(N) a 

(di,d 3 )=l (B,di«fe)=l 

The character sum S^x, 3, /i, <^i, ^2) is multiplicative in x, 9, ^ in the following sense. If g = <7i<72 and 
/i = h\h 2 are any decompositions such that (gihi, g 2 h 2 ) = 1, and correspondingly x = X1X2 where 
Xi is a Dirichlet character modulo gihi, then 

S(x> g, h, d ll d 2 ) = S(xi,gi, h 1 ,d 1 ,d 2 )S(x2,92, h 2 , d l ,d 2 ). 

Indeed, this follows easily upon writing the summation variable d mod gh as 

d = diei+d 2 e 2: dimodgihi, 

where e\ and e 2 are fixed integers such that 

ex = 1 (mod gihi), e% = (mod g 2 h 2 ); 

e 2 = (modgihi), e 2 = 1 (mod g 2 h 2 ) . 

In other words, we have a decomposition over the primes 

S(x,9,h,di,d 2 ) = Y[S(x P ,g P ,h p> di,d 2 ), 
p 

where the subscript p denotes the p-part. 

Let us fix p for a moment and use the notation 

g p =p J , h p =p s , N p =p u . 

Note that Xp is a primitive Dirichlet character modulo p v . 

Case 1. If v = (i.e. p\ N), then x P is trivial, so that 

fl, 6= 0, 

S{Xp,9p,h p ,di,d 2 ) = E 1 = \ P 5 ~P & ~ l i s > and p | did 2 p 7 , 

d(p T+,s ) Ij?*-2j7 s - 1 , <5 > andp|did 2 p 7 . 

(did±d 2 ,p 7+,5 )=p T 
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For 5 = the right hand side follows by observing that there is a unique d satisfying the condition, 
since 7 > implies p \ d\d2- For 5 > and p | d\d2 the right hand side follows by observing that 
7 = and exactly one of d\ and di is divisible by p, hence the condition on d is automatically 
satisfied. For 8 > and 7 > the right hand side follows by observing that p { g?ic?2 , so that the cf s 
satisfying the condition are in bijection with the reduced residues modulo p s . For S > 0, 7 = 0, and 
P \ d\d2 the right hand side follows by observing that the condition on d is (d\d± c?2,p) = 1 and 
there are precisely p s_1 reduced residue classes modulo p 5 that do not have this property. 

Case 2. If v > (i.e. p | TV), then p \ d±d2, 7 + 5 ^ v , and \p induces a nontrivial character modulo 
pi+s Q f conductor p v ', so that 

X P (Tl)S(x P ,g P ,hp,di,d 2 ) = E X P {Td2did)= E Xp {d) 

d(p-<+ 5 ) d(p T+5 ) 
(did±d 2 ,p T+5 )=p T (d-l,p"'+ a )=p^ 

1, (5 = and 7 ^ v, 

p s -p 5 ~ 1 , 8 > and 7 ^ j/, 

— p* 5-1 , (5 > and 7 = av — 1, 

0, S > and 7 < z> — 1. 

For <5 = and 7 ^ ^ the right hand side follows by observing that there is a unique d satisfying 
the condition in the second sum, and this d is congruent to 1 modulo p v . For S > and 7 f the 
right hand side follows by observing that the d's satisfying the condition in the second sum are in 
bijection with the reduced residues modulo p s , and all these d's are congruent to 1 modulo p u . For 
S > and 7 = ^ — 1 the right hand side follows by observing that the condition on d only depends 
on d modp u , hence in this case 

E* xM)=v 5 - 1 E* %>(<*), 

d(p~<+ 5 ) d(p") 
(d-l,p~< + s )=p~i (d~l,p v )=p u - 1 

and the sum on the right hand side equals — 1. Indeed, for v = 1 this is obvious, while for f > 1 it 
follows from the fact that u 1— > u + 1) is a nontrivial additive character modulo p. Finally, 

for S > and 7 < f — 1 the right hand side follows by observing that the second sum does not 
change when we multiply it with a complex number of the form X P (p v u + 1) 7^ 1 ■ 

Our findings imply that 

S{Xp,9p,h p ,d 1 ,d 2 ) = S(x P ,9p,hp,di,p,d2 lP ), 
hence we have an Euler product decomposition 

p 

where for p \ N 

E r,p( s ) : = T _ l-X-p E ^STflA E 

y a, 0^0 F 7,<5^0 1 

min(a,/3)— min(7,a+/3)— 

while for p | AT 

3^( S ):=X P (Tl) E ^W 5 (WV,M). 

7,<5>0 ^ 

Now we insert the above calculated values of S (x P > P 1 , P 5 > P a > P^) ■ 
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Case 1. If p \ N, then we obtain 

oo ^ 2 °° 1 -^oooo _^ 

(1-P ' )-^ p (s) = ^ — + — — 5(A-1) + X] X] d7m +5(A-1) 

7=0 ^ ^ 5=1^ ^ 7=U=1 

Indeed, the first line contains the contribution of a = /3 = 0, and the second line contains the rest 
(where 7 must be zero). We sum all the geometric series: 

(i- P -*-»)zU s ) 

p- 2 , p-i , 2 ^ , p-i A 



l-p-M p A (l-p 1 - A ) p A +^(l -p-M) p A (l-p~ A ) \ p A (l-p 1_A ). 

(1 -p 1 -*) +p~ x (p- 2)(1 -p~ M ) + p~ A ^(p- 1) 2p" A 

(1 -p 1 - A )(l -p-^) + 1 -p^ A 

1 - 2p~ A + p- x -^ 2p- x - 2p- x -^ 



(1 - _p 1 - A )(l - p-^) (1 - p^ A )(l - p-») 
1 -p- x -f+ 



(1 -p 1 - A )(l -p-^ 1 ) ' 
Hence for p j AT we have 



;=:, r,p( s ) — 



(1 -p 1 - A )(l -p-A 1 )' 

Case 2. lip\N, then we obtain for p u \\ N 

Xp(Tl)" r , p (s) = 2J " i+(„_i) M X p 5(A-l) ^ X X jj7/*+<5(A-1) ■ 

7=v 5=1 7=1/ J=l 

We sum all the geometric series: 

Xp(,Tl)-r,p( s ) - -p-M) ~ p A +("-l)M(l -p!- A ) + p A +"^(l -p!- A )(l -p-^) 

_ (1 -p 1 ^) -p~ A +^(l -p"^) +p~ X (p - 1) 
p^(l -p 1 - A )(l -p~») 
1 -p- A +^ 



P^(l -p!- A )(l -p-' 1 )' 

Hence for p | N we have 



1 - p- A +^ 



(1 -p!- A )(l -p-^)' 

Collecting the above results we arrive at (61). □ 
4.3. Applying Voronoi summation. We substitute (59) back into (54) getting 



^ (S,N) = 1 
powers of two (d 1 ,d 2 ) = l 



Mi,M 2 ,G>l ± (S,N) = 1 N\c Yy L Z ' 



(62) 

V^' (9j+ a '2\ c Afi,Af 2 /-r 7 \ 



gh=c f(dld 2 h) 
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where 



S d u d 2 ,c(f>9,h) := 2^ A s (mi)Ag(m 2 )e I ±mi — + m 2 — I -F dl) j 2;C (mi, m 2 ), 

"11,11! \ / 

{X > y) = (xyW 1 " 3 [ Tl0& dp 8^5*' J e [ T d^-c 

By applying Proposition 2 for the summation variables mi and m 2 , we see that 

dxd 2 h 2 S^^(f,g,h) 
is a sum of terms (suppressing M\ and M2 from the notation for simplicity) 



(63) 



with 



X g (m 1 )X g (m 2 )e (^ mi -^-^ =f (m 2 fmij^j ^dt',d 2 ,e( TO i> TO 2) 



(64) F dud2 Jm um2 ):=j o ^ F di ^ c (*, y) J g [~^~ ) J g {^J^ ) ** *y. 
If g = EV is an Eiscnstcin series, there are three additional polar terms 



(65) 
(66) 



A(m 2 )e (rm 2 ^H ^ ^JiJa.cC^ TO 2) 



/"CO /*00 

(67) / / F^c 2 (x,y)P^ 2h (x)P r % lh (y)dxdy, 

Jo Jo 

with 



poo 

(68) ^t:*,c("»i.y)== I d^)^ \ 

(69) m 2 ) := [ (*, „) J* (^^) 



(70) f£(t) 



C(l±2ir)(t/c 2 ) ±8r , forr^O, 
log(i/c 2 ) + 27, for r = 0. 



We proceed to analyze the four terms (63), (65), (66), (67). It will turn out that the first 3 terms 
are small, but (67) contributes to the main term. 

4.4. The contribution of (63). We observe first that in (64) the arguments of the Bessel functions 
are large: 

-2 7^2 7^2 



> -V^r > ^ . ~ —. 7TTT, 7T77 > 1=7 



d 2 h d 2 G C s d 2 y/M 1 M 2 C s (d 1 d 2 ) 1 l 2 {M 1 M 2 ) 1 / 4 C £ (T + K) 

by (55), (57) and the fact that h ^ c. A similar estimate holds for Air^m 2 y/ {d\h). In view of (9) 
this is large. In particular, by the rapid decay of the Bessel X-function (32), among c (mi,m 2 ) 
we only need to consider F£'~^ c ( m ij m 2) as the contribution of the other 3 expressions is negligible 
(or zero). 

We show now that the sum (63) can be truncated efficiently. Fix any y € [(1/2)M 2 , (5/2)M 2 ] in 
the integral defining F di '~j^ 2 c (mi, m 2 ). By (56), we can restrict the x-integration to 

(71) d 2 1 x = 4y{l + 0(T^ 1 )) 
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at the cost of a negligible error. In this range we have 

d j . . ~ s ( T di 

jF du d 2 ,c(x,y) <j C 



fai'm'^-"" ~ y Ml Td 2 c 

so that by Lemma 3 the integral is negligible unless 

T di \ \ZMld 2 h 



> c- 



Mi Td 2 c ) \lfn\ 



that is 

(72, mi ,M,.:=^(M! + ^ 3 

Here we used (57) and the fact that gh = c. Similarly, we can assume 



(73) to 2 < M 2 * := C e d\ 
For convenience we observe, by (57), 

(74) M* := max(Af*, M*) x C e d x d 2 
To summarize, 



Mi 



( iThf 

\ M 2 ' (Tg) 2 

(Th) 2 (M^h) 1 ' 2 



{M X M 2 )V* (Tff) 2 



<^(/>5» = ^ E A ^) e (j£g) E Vm 2 )e(Wm 2 )^- 



(75) 



oo /-oo 



^0 



up to negligible error terms and the contribution of the three polar terms (65)-(67) that we discuss 
in a moment. 

If we substitute this back into (62), then the summation over / produces the exponential sum 

f(dld 2 h) 
(g^fd,2 f ,d\ d,2 h)—d\ 

We estimate the double integral in (75) using (55), (71), and (31) as 

< (TMiM 2 ) mm , -rjj < C — -rjj. 

\ d\ d\ ) (m,im 2 ) 1 ' 4 T(toito 2 ) 1/4 

By Lemma 6 and the support of uj, the contribution of (63) to (62) is 

«TKC £ V 1 Qh (™i ~ ™-2)did 2 ) |A g (mi)A g (m 2 )| 

«TKC e V 1 V- 1 V/ V l A g( TO i) A g( m 2)l 

^(ditfe) 3 / 2 (gh)V 2 hT V ^ . (mim 2 ) 1/4 

£|(mi- m 2 )did2 

By (33) the innermost sum is 

< E 1 f i A ^f + ^f ) « (m-)v^ f i + 

^ \ 771/ 771, / \ " 

t\{mi —m2)did 2 
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hence in the end the contribution of (63) to (62) is, using also (55) and (74), 
(76) 



<r TKC £ V 1 f g 3/2 | G 1/2 (MiM 2 )V 4 T 2 ^/ 2 ( Ml M 2 ) 3 / 4 \ 

4- ^ Udid2)(MiM 2 )V4 + ( dl d 2 )T 2 + (M x M 2 )3/4 + T 4 ; 

<r TKC e V 1 f ( M * M 2) 1/2 . W(! , T 2 {M\M 2 ) 1 ^ 2 (MiM 2 ) 3 / 4 \ 
^ 6 2 I (d^K 3 (dida)T*K K 5 T* J 

1 /(7 1 / 4 c 1 / 4 t 2 ^ 1 / 4 C 3 / 8 \ 

< Tif C e > „ , , — r + 777777 H h -7777- 

<5 2 d 1 d 2 sCC 1 /4+£ \ / 

, .it / T 4 A 3 

4.5. The contribution of (65) and (66). We show that the integrals (68)-(69) are negligible in 
the ranges (72)-(73). Starting with the first, we see by the rapid decay of fl 3 in the definition 
of F dl ^ d2 ^ c (x,y) that the contribution of \d\x — d\y\ ^ T e ~ 1 d\Mi is negligible. Introducing w := 
(y/ dfx — yj d\y) 2 and applying several integrations by parts with respect to w shows that the 
contribution of yj wd\x did 2 cT 1+E is also negligible. Using \d\x — d\y\ <C \Jwd\x we infer that 
we can restrict the integration to 

\d\x-d\y\ ^Z:= & min (^^-,did 2 cT 

at the cost of a negligible error. In other words, writing z := d\x — d\y we can approximate the 
above integral by 

1 f z [°° ( d\x-z \ ± (djx-z\ , ( 4^/rm \ 

with negligible error. For \z\ ^ Z we have 

ft* \**'*» C p d 2 J ^'dih ^ d 2 ) } «J ° ^ M 2 + Ml + ( d 2 M 2 )(did2c) j 

whence by Lemma 3 and (55) the last integral is 

^ cz \7m) <<j z {(m 1 m 2 )^j << ^ cz { w 

/cVsV _ fT + K 



K 2 V A 2 



Choosing j g N sufficiently large and using (9), (55), (72) we conclude that the polar term (65) is 
negligible. In the same way we see that the polar term (66) is negligible. 



4.6. The final polar term. It remains to estimate the contribution of the polar term (67). This 
requires some non-trivial manipulation. First we remove the dyadic decompositions by summing 
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over M\,M<z,G using (53), but we keep in mind that the decay of O3 allows us to restrict to 
<5, d\, c?2, c ^ K £ , up to a negligible error. This gives 

(77) 



2 ^ <5 2 did 2 V <£(did 2 c) cV2 
± (<5,A r ) = l iV|c V ; gh=c 

(di,<fe)=l 



f(dfd 2 h) 



d\d 2 h 2 \ d\ 



oo poo 



^0 



^3 Tlog 



d|y ±K 2 c 5 A d 2 d 2 x y \ ^ ( y/djx - y/d$jj) 
d 2 x' 8ny/xy' J~£< 



did%c 



E, . . 7->i / \ dx dy 
F r,d 2 h\ x )K,dih\y) ( xy )3/4 ■ 



up to a negligible error coming from the fact that X)j>o w ( a; /2 : ') is n °t necessarily 1 for x < 1. In 
order to simplify the notation, we will only consider the (+, +)-tcrm in the last sum and drop the 
superscripts at P. 

We observe that the <5, d\, d2-sum is rapidly converging due to the decay properties of O3. A trivial 
estimation shows that the double integral is 0(C £ ). Now we replace O3 by O3 which introduces by 
(52) an admissible error of 

(78) 0((TK) 1+£ (T + K)- 1 ) . 

Next we make a change of variables 



z := xy, 



4v , 

w := -~ 1. 



dfx 



By Taylor's formula, 



dfx 



Tlog-|^ = Tu; + 0(Tu; 2 ) =Tw; + 0(T £ ~ 1 ) 



in the range where is not negligible. Again by Taylor's formula 

2 



hence 



(Jd^-^/dfy^j =^d 1 d 2 Vzw 2 (l + 0(w)), 



d\d 2 c 



4c 



(1 + ocr- 1 )) 



in the range where O3 is not negligible. Therefore the double integral equals 



00 />oo 



1 JO 



3 \^Tw 

x Pr,d 2 h 



±K 2 c 8 i d\d 2 z\ e f y/zv? 



d 2 y/z 



<C ) 
Pr,d x h 



4c 



di y z(l + w) \ dw; dz 



(79) 



00 />oo 



00 JO 



3 Tw 



di V 1 + w 
±K 2 c S 4 d 2 d 2 z 

8lTy/~Z ' yJ~Q 



2(l + w)z 3 / 4 



C^T 6 " 1 ) 



XPr,d 2 



d 2 yfz 
di 



P, 



r,d\h 



d\Jz 



dw 



dz 

2z 3 / 4 



+ 0(T £ - i ) 



We rewrite the last w-integral using the definition (51): 



(80) 



00 p 00 



Oo x 



±K 2 c S 4 d 2 d 2 2 



exp iTwx =F 



~2c 



dx dw. 



AN ASYMPTOTIC FORMULA 



25 



This double integral is not absolutely convergent, so we consider 



oo JO 



±K 2 c S 4 d 2 d 2 z \ / e 



(81) L .= / f^2 x, == — exp iTwx dxdw. 

87r\/i' / V 2c 



for e > 0. The x-intcgral in (80)~(81) is a Fourier transform decaying rapidly in w, therefore by 
cutting the w-integral at larger and larger parameters we see that limsup e _ > .Q + \I — I E \ is smaller 
than any positive number, i.e. I = lim e _>o+ Ie- As the double integral (81) is absolutely convergent, 
we can change the order of integration there and compute the w-integral using [GR, 3.323.2]: 

Here the integrand is rapidly decaying in x, hence by a limsup argument as before we see that 
7= to U -f<* f^)"* ffiJ^^.^lexpf^U. 



mi J £ i i i u sj^ I iAj . , . I v^./vjj I 

>0+ \V Z J Jo \ S^y/Z \l (j I \ 27T\/2 

Using also the definition (49) we can summarize that the w-integral in (79) equals 



We integrate this over z and substitute it back into (77) getting 



V ^ (Wid 2 ) 2 ^ 0(did 2 c) ^ /i 2 ^ X V di 

± (5,iV)=l V 21 N\c^ y ' gh=c f{d\d 2 h) V 

(dud 2 ) = l (g,did 2 ) = l ( g f T d 2 ,dld 2 h)=di 



00 r°° r°°- ( S 4 d 2 d 2 z\ / ic ( „ 2 y 2 K 2 
fii ( y, LJ- ] exp ( ( ^ 2 T 2 + ^- 



o Jo Jo \ \l C 



"r.dih — -, — \dydx- 



di J ' 1 V. ^2 ; 2z 

Let us consider the case r ^ and insert the + case of (70). Then we can recast the preceding 
display as 

(82) 



± (<S,jV)=l V ^ JV|c ^ V ; 9 h=c f(dld 2 h) V 

Ox *, y, exp f ± ( ^ T + V ) ) ( ^ j ^ 



'o Jo Jo \ y C J 

We make a change of variables 

c 

— 7= = v > 



and write the triple integral as 

?4 J2 j2„2 \ / 1 \ / „ \ 2 "" 



fOO fQO />oo 



/°° x, y, S -^= exp f±ii« ((xT) 2 + (i^) 2 )) (—^j-rs) dydx 
/o Jo Jo \ n 2 v 2 VC J V 2 ) Knvd^h 2 J v 
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Finally we insert the definition (48) and arrive at 



oo />oo poo 



JO JO 



Wl{x)W2{y) ± [ Gq(s) '"t^T ^ 1 ds 



x exp (±- W ( (a r)» + (I^T) ) ) dyd X - 

by Mellin inversion. We compute the w-integral by [GR, 3.381.4]; the change of the order of integra- 
tion can be justified similarly as before. Reorganizing, we obtain 

)2ir .oo poo 
l I W l{x) W.{y){{ X T? + { \yK?)- dydx 

1 ( N 



x ™ y (2) {^dVc ) T{2s - 2 ^ cxp ^^ - ir » ds ' 

where by (8) the x, y- integral is simply 

We substitute this into (82), recall the definition (60) and recast (82) as 

C(l + 2ir) 2 TKM ir (T, K)— G (s) EJs) N 2s ~ 2lr 2(2n)- 2s+2lr T(2s - 2ir) cos(7r(s - ir)) ds, 

2m ./(2) 



where ^. r (s) = S^(s) denotes the function in (61). 

To summarize, the contribution of (67) equals, up to an admissible error, 

(83) ((l + 2ir) 2 l[(l-p- 2 - 4lr )TKM lr (T,K)^- j G (s)Z r (s)ds 

p\N m ^ (2) 

with the kernel 

Z r (s) := 2(2tt)- 2s+2 "T(2s - 2ir) cos(tt(s - ir))((2s - 2ir)((l + 2s + 2ir). 
Using [GR, 8.334.2 & 8.335.1] the kernel equals 

Z r (s) = ^- 2s+2 * r J^"" -M 2s - 2ir)C(l + 2s + 2ir), 
T(i - s + ir) 

hence by the functional equation for the Riemann zeta function we obtain 

Z r {s) = C(l - 2s + 2ir)C(l + 2s + 2ir). 

Let us assume r ^ 0, then the integrand in (83) is an odd function of s which is holomorphic 
except for a simple pole at s = as well as possible poles at s — ±ir. It follows that the s-integral 
equals half the sum of its residues, that is, (83) equals 



(84) 



i«l + 2w) 4 [](1 p- 2 - 4 "-) TKMir(T, K) 

p\N 

-|C(1 + 2zr) 2 C(l + 4^r) [](1 - p- 2 - 4lr ) G {ir) TKM lr {T, K). 



p\N 



Here we also used that Go(— ir) = — Go(zr). The (— , — ) case in the last sum of (77) gives a similar 
contribution to Ai-i r {T, K ), so that the second line of (84) together with the corresponding part of 
the (— , — ) case precisely cancels (46). The (+, — ) and (— , +) cases contribute to £o- 

The case r = can be treated similarly and gives a linear combination of £qj £ij £2 in Theorem 3. 
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This completes the discussion of the off-diagonal term. The various error terms (50), (76), (78), 
(79) encountered so far are admissible for Theorems 1-3, while the first line of (84) and its counter- 
part with ir replaced by — ir have the desired shape for Theorem 2. The proofs are complete. 

5. Appendix 

In this Appendix we deduce Lemma 1 from the modularity of g = E r (r ^ 0), in order to 
emphasize the analogy with the cuspidal case. We follow closely [HM, Section 2.4]. 
By [Iw3, (3.29)] we have the Fourier decomposition 

6 (| + ir) g(x + iy) =6{\ + ir) yi +ir + 6 (§ - ir) y^"' 

oo 

+ Xg(n)Ki r (2wny) cos(27mx), 

n=l 

where 

(85) 9(z) tt- z T(z)C(2z) and \ g (n) := ^ (| 



For convenience we introduce 



D ±1 {g,x,s) := ^D(g,x,s) ± ^D(g,-x,s), 



i.e. 



D +1 (g, x, s) = A 9 (n) cos(27rna;)n 

n=l 
oo 

D^ 1 {g 1 x 1 s) — ^""^ \ g (n)i sin(27rnx)T 



n=l 

It will be more pleasant to work with the Maass shift (cf. [DFI, (4.3)]) 



( d \ 

g(x + iy) := -y [i — + —J (\ + ir) g(x + iy) 



which is a weight 2 Eisenstein series of Laplacian eigenvalue 1/4 + r 2 . By [DFI, Section 4] we have 
the Fourier decomposition 

g(x + iy) = g const (y) + g SCI (x + iy), 

where 

(86) .g const (y) :=-(£ + ir) 9 (± + ir) y? +ir - (± - ir) 6 (± - ir) yi~ ir , 

(87) g^x + iy)- - X9{n) 



■= X] 9 r {^ir^ny) coa(2nnx) + V 2 l r {^ny)i sin(27rnx)} 
, v ri 



and V±l is as in [DFI, (8.27)] 



2,ir 

v 2 %(y) ■■= WiMW T Q + r 2 ) W-x^y). 

Utilizing the functional equation 



~ / a iy\ _ (a i \ 

<?- + -= -<7 — + — > y > o, 

\C c J \ c cyj 

which is a consequence of modularity, we sec by standard estimates that 

(89) g scr + < 9 , c min^-^^e- 2 ^). 
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Therefore, taking Mellin transforms we obtain by (87) 

(90) J™ §~ + f ) f = ^ (3 "F (g, ±a) , Ms > 1, 
where 

(91) F(g, x, s) := S+ 1 (s, i r).D+ 1 ( 5 , x, s) + $ 2 " 1 (a, Z r)^" 1 (g, z, s) 
and 



Jo 2/ 



The last function is the same as [DFI, (8.25)] except that we deleted the 4 from the denominator. 
The reason is that we would like to apply [DFI, Lemma 8.2] but that lemma requires this correction, 
because the initial identity [DFI, (8.30)] is missing a factor 1/4 on the right hand side. By [DFI, 
Lemma 8.2], 

(92) ^ {s> ir)=( s -i)r( S -±^)T 



2 

(93) $2 -WH2r(i±f±i)r^-- + 1 

We shall derive the analytic properties of D(g, x, s) from those of F(g, x, s). Splitting the integral 
in (90) and applying (88) in the form 

( a iy\ {® *\ ( 1 \ (V\ 

5ser - H = -ffser 1 ~ Sconst — ~ ffconst ( ~ , 1 > J/ > 0, 

\c c ) \ c cyj \cyj \cJ 

(90) becomes 

(94) £)'f (g, " s) =V~c f (&- + y*- 1 * - 5- (-? + ^ ^ + ^ C «), 
where 



7T/\C/ Jil\ cc / \ cc / J !/ 



(95) P(ir,c,s) := - Vc j<?const (~) + Sconst Q 



y 2 — ■ 

y 



By (89) the first integral in (94) is absolutely convergent for any s € C, hence it defines an entire 
function. Moreover, it becomes its own negative under the substitution - — > — - and s — > 1 — s : we 
shall say it is symmetric for short. The second integral (95) can be calculated explicitly using (86), 

(96) P(ir, c, s) = £ (i ± ir) (| ± ir) - j^-) , «a > 1- 

The function P(ir,c,s) is meromorphic and symmetric in s £ C, hence (^) S F (g, -, s) is also 
symmetric and differs from P(ir,c,s) by an entire function. If we apply this conclusion for — - in 
place of | and combine (91) with D ±1 (g 1 — x, s) = ±D ±1 (g,x,s), then we see that the following 
functions are entire and symmetric: 

(97) (~) S $t{sMD +1 (g,- c ,s)-P{ir,c,s) and ^\s, ^D" 1 (g, ~, a 

Now we can argue exactly as on [HM, p. 597] (with k = 2) to see that D(g,x,s) satisfies the 
functional equation (29), cf. [HM, (15)]. Here we note that the functions ^^ it on [HM, p. 597] 
should have been halved, therefore in [HM, (19)-(20)] the factors 2 2s are really 2 2s_1 . 
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Wc also need to analyze the poles of D (<?,§, s) for which we go back to (97). We see by (93) 
that -D -1 (g, f , s) is entire. By (92) and (96) the poles of D +1 (g, ^, s) are at s = 1 ± ir; they are 
simple with residues (cf. (85)) 



(- 

V c 



$2 (1 ± ir, ir) 

We note that the potential poles at s — ±ir coming from the poles of P(ir, c, s) are canceled by the 
poles of $2 _1 (s,ir) there, while the potential pole at s = 1/2 coming from the zero of "J^^s, zr) is 
canceled by the zero of P(ir, c, s) there. This concludes our proof of Lemma 1. 
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